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Abstract 

We investigate the production of pions in heavy-ion colhsions in the energy range 
of 1 - 2 GeV/A. The dynamics of the nucleus-nucleus collisions is described by a set 
of coupled transport equations of the Boltzmann-Uehling-Uhlenbeck type for baryons 
and mesons. Besides the A^(938) and the A(1232) we also take into account nucleon 
resonances up to masses of 1.9 GeV/c^ as well as tt-, rj- and p-mesons. We study in 
detail the influence of the higher baryonic resonances and the 27r-production channels 
{NN — > NNtttt) on the pion spectra in comparison to tt^ data from Ar + KCl colli- 
sions at 1.8 GeV/A and 7r°-data for Au + Au at 1.0 GeV/A. We, furthermore, present 
a detailed comparison of differential pion angular distributions with the BEVALAC 
data for Ar + KCl at 1.8 GeV/A. The general agreement obtained indicates that the 
overall reactions dynamics is well described by our novel transport approach. 



1 Introduction 

Relativistic heavy-ion collisions provide a unique possibility to study nuclear matter at 
high density and far away from equilibrium. During the course of a heavy-ion collision at 
1-2 GeV/A the nuclear matter is compressed up to 2 — 3 times normal nuclear matter 
density pQ before it expands again. The experimental probes that provide information 
about the compressed stage of the reaction are collective observables such as flow patterns 
or differential spectra of particles that have been additionally produced during the collision. 

Over the last 10 years transport theories as BUU ||, |, |] and QMD, IQMD @, |, |, § 
have been very successful in describing the reaction dynamics of heavy-ion collisions. It has 
been found that the experimental meson spectra can be well understood when assuming the 
excitation and subsequent decay of nucleon resonances during the compressed stage. Since 
pions couple strongly to these resonances, the differential pion spectra provide a well suited 
probe for the dynamics of baryonic resonances. Furthermore, due to the low production 
threshold, pions are produced and reabsorbed quite frequently and thus provide a signal 
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for the whole dynamical evolution of the heavy-ion reaction. This implies in particular 
that a transport theoretical description of heavy-ion collisions has to reproduce the pion 
yields correctly before one can draw any further conclusions on more specific channels from 
the model. Especially for observables like dileptons the pion annihilation (7r"'"7r~ — > e"'"e~) 
plays an important role; also the production of kaons (via vrA^ KA) and ?7-mesons (via 
ttN — > rjN) is strongly influenced by the pion induced channels. 

In this work we present a new implementation of the BUU-model, denoted by 
Coupled-Channel-BUU (CBUU), which as an extension to previous realisations in ^ - 
or compared to the transport model IQMD ^, - includes all baryonic resonances up 
to masses of 1950 MeV/c^. The higher resonances are expected to contribute to the low 
as well as to the high momentum regimes of pion spectra via Ivr and 27r decay channels, 
respectively. Apart from the pions (vr"*", vr*^, vr") we now also explicitly propagate ?7's and 
p- mesons {p^, p^, p~) as mesonic degrees of freedom thus including hadronic excitations 
up to about 1 GeV of excitation energy. After introducing our model in section 2 and 
the various elastic and inelastic cross sections in section 3 we present detailed compar- 
isons with 7r~data taken at the BEVALAC |1C] and the SIS in section 4. A summary and 
discussion of open problems concludes the paper in section 5. 



2 The CBUU-Model 

2.1 Basic Equations 

In line with refs. [Q, |2|, ^, 11] the dynamical evolution of heavy- ion collisions or hadron- 
nucleus reactions below the pion-production threshold is described by a transport equation 
for the nucleon one-body phase-space distribution function fN{f,p,t), 



dt 



+ |— + VpU{r,p)\Vrf{r,p,t) 

+ {-^^%^ V,t/(f, V,/(f, p, t) = hoiAfir, p, t)], (1) 



where r and p denote the spatial and the momentum coordinate of the nucleon, respec- 
tively, while N stands for a proton (p) or neutron (n). The effective mass m*{f, p) in 
eq. (|^) includes the nucleon restmass rn-jy (= 938 MeV/c) as well as a scalar momentum- 
dependent mean field potential [/(f, p), 

m*{f, p) = rriN + U{r, p). (2) 

The nucleon quasi-particle properties then read as 



E = ^m*{f, p)'^ +p^. (3) 

Physically the l.h.s. of eq. (|l]) represents the Vlasov equation for a gas of non-interacting 
nucleons moving in the scalar momentum-dependent mean- field potential U{'r, p). The 
Vlasov equation as given above can be derived from a manifest covariant transport equa- 
tion using scalar and vector self-energies depending on the four- momentum (p^^ = {E, p)) 
of the particles [11, when neglecting the vector self-energies. Due to the quasiparti- 



cle mass-shell constraint given by (y) the scalar potential effectively depends only on the 
three-momentum of the nucleons as {U{p, E{p)). 

The r.h.s. of the BUU-equation (i.e. the collision integral Icoii{f{r, p, t)) ||, |2|) describes 
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the time evolution of /Ar(f, p, t) due to two-body collisions among the nucleons. For ex- 
ample, the alteration in the one-body phase-space distribution function f{fi, pi, t) due to 
the elastic scattering of two nucleons 

► N3 + iV4, 



Ni+N2^ 
with momenta {pi, i = 1, ..,4) is given by 



Icoll [fi{r, pi,t)] 



E 

2,3,4 
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(fp2 J d^P3 J dn^vu 



X ifsir, P3, t) fi{r, Pi, t) /i(f, pi, t) /2(f, p2, t)) 

- fi{r, Pi, t) /2(f, p2, t) fsif, p3, t) fi{r, p4, t)), (4) 

where '^'^^^^^'^ is the in-medium differential nucleon-nucleon cross section, fi = l — fi{i = 
1, .., 4) the Pauli-blo eking factors and V12 the relative velocity between the nucleons A^i and 
in their center-of-mass system. The factor 5 = 2 in (Q) stands for the spin degeneracy 
of the nucleons whereas 3 4 stands for the sum over the isospin degrees of freedom of 
particles N2,N3 and A^4. 



For energies above the pion-production threshold one also has to account for inelastic pro- 
cesses such as direct meson production channels or the excitation/deexcitation of higher 
baryon resonances. In the CBUU-model - to be described here - we explicitly propagate 
the mesonic degrees of freedom n, r], p and a scalar meson a that simulates correlated 
27r pairs in the isospin 0-channel. Besides the nucleon and the A(1232) we, furthermore, 
include all baryonic resonances up to a mass of 1950 MeV/c^: i.e. iV(1440), iV(1520), 
A^(1535), A(1600), A(1620), iV(1650), A(1675), iV(1680), A(1700), iV(1720), A(1905), 
A(1910) and A(1950), where the resonance properties are adopted from the PDG |13U . 
Introducing one-body phase-space distribution functions for each particle type leads to 
equations similar to eq. (|^) for each hadron. Since the different particle species mutually 
interact the integro-differential equations are coupled by the collision integrals. Schemat- 
ically one can write down the set of coupled equations in the following way: 
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(5) 
(6) 

where Df abbreviates the left side of the Vlasov-equation. The collision integrals on the 
r.h.s. of eq. (P) have formally the same structure as the one given in eq. (^). In addition to 
elastic scattering processes they now also contain all allowed transition rates. Denoting the 
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nucleoli by N and the baryon resonances listed above by R and R' , we include explicitly 
the following channels: 



• elastic baryon-baryon collisions 

NN ^ NN 
NR ^ NR 

The elastic A^A^-cross section is described within the parametrization by Cugnon 
|[l|, while the cross section for elastic A^i?-scattering is evaluated by using an 
invariant matrix element which is extracted from the A^A^-cross section (cf. Section 
3). For A^i?-scattering we also allow for a change in the resonance mass according 
to the corresponding distribution function (cf. Section 3). 

• inelastic baryon-baryon collisions 

NN ^ NR 
NR ^ NR' 
NN ^ A(1232)A(1232) 

For the cross section of the reaction A^A^ A^A(1232) we use the result of the OBE- 
model calculation by Dimitriev and Sushkov |jl5|. In order to obtain the A^A^ <-> 
NR-cioss sections for any of the higher resonances we exploit the resonance model 
described in more detail in Section 3, while the parametrization for the A^A^ <-> 
A(1232)A(1232)-cross section is adopted from Huber and Aichelin 

• inelastic baryon-meson collisions 

R ^ Ntt 
R A^vrvr 

^ A(1232)7r, A^(1440)7r, Np, Na 
Af(1535) ^ Nr] 
NN ^ NNtt 

Besides the production and absorption in baryon-baryon collisions, baryonic reso- 
nances can also be populated in baryon-meson collisions and subsequently decay to 
baryons and mesons again. A detailed discussion of the respective cross sections 
and decay widths is given in Section 3. Here we note that the 27r-decay of higher 
resonances is modelled via subsequent two-body decays as indicated above. 

• meson-meson collisions 

p TTvr (p-wave) 
a vrvr (s-wave) 

For the pure mesonic cross sections we use the Breit-Wigner parametrizations given 
in Section 3. 
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2.2 The Test-Particle Method 



The CBUU-equations are solved by means of the test-particle method, where the 
phase-space distribution function /Ar(r, p, t) (e.g. for nucleons) is represented by a sum 
over (^-functions: 

, NxA 

fir, P,t) = - ^^(t)) X ^ (P- • (7) 

i=l 

Here denotes the number of test-particles per nucleon while A is the total number of 
nucleons participating in the reaction. Inserting the ansatz (^) into the CBUU-equations 
dll) leads to the following equations of motion for the test-particles: 

dt ~ E ^'P^i^^'P'Wj 

dpi(t) m* -± ^ , s.. 



dt E 

Again these equations of motion are consistent with those from a covariant transport 



equation |11| when using only scalar self-energies for the baryons. Thus the solution of the 
CBUU-equations within the test-particle method reduces to the time evolution of a system 
of classical point particles according to eq. (^). For the actual numerical simulation we 
discretize the time t and integrate the equations of motion employing a predictor-corrector 
method |17|. We want to note that in our model pions are treated as 'free' particles, except 



for the Coulomb interaction. It remains to be seen if pion selfenergies as suggested in refs. 
PI will alter our results. 



2.3 The Colhsion Integrals 

The collision integrals occurring in eq. (P) contain either particle-particle collisions or the 
decay of baryonic or mesonic resonances. 

For particle-particle collisions we employ the following prescription: The test-particles 
collide with each other as in conventional cascade simulations with reaction probabilities 
that are calculated on the basis of free cross sections. The numerical implementation 
additionally accounts for the Pauli-blo eking of the final states while the collision sequence 
is calculated within the Kodama algorithm [^] (cf. ref. |^]) which is an approximately 
covariant prescription. 

The decay of a resonance, furthermore, is determined by its width T{M). However, during 
the course of a heavy-ion collision the resonance may also decay due to collisions with other 
particles (eq. NR — > NN). This coUisional broadening of a resonance is described by the 
collision integral for particle-particle collisions. Here we will concentrate on the numerical 
method used to account for the first decay mechanism. All resonances treated in the 
CBUU-model are allowed to decay into a two-particle final state, i.e. in every timestep 
of the simulation we calculate the decay probability P for each resonance assuming an 
exponential decay law 

P = l- e-r(M)/h^At^ 

where At is the time step size of our calculation, T{M) is the energy-dependent width 
of the resonance and 7 the Lorentz factor related to the velocity of the resonance with 
respect to the calculational frame. We then decide by means of a Monte Carlo algorithm if 
the resonance may decay in the actual timestep and to which final state it may go. If the 
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chosen final state contains a nucleon, which e.g. is Pauli-blocked, we reject the resonance 
decay. 



3 Explicit numerical Implementations 

3.1 Nuclear Mean-Field Potentials 

Particles propagating inside nuclear matter are exposed to the mean-field potential gen- 
erated by all other particles. From Dirac-phenomenological optical-model calculations 
|11, 21| it is known that elastic nucleon-nucleus scattering data can only be described 



when using proper momentum-dependent potentials. Here we employ the momentum- 



dependent mean-field potential proposed by Welke et al. |22], i.e. 



U-if, p)=aJ^ + B^^ + 2^ / d^'^p^. (10) 
Po Po Po J 1+1 p-p"' 



A 

As an extension of the momentum-independent Skyrme type potentials for nuclear matter 
1^ ^ the parametrization ([lO|) has no manifest Lorentz-properties. However, definite 
Lorentz-properties are required for a transport model at relativistic energies. To achieve 
this goal we evaluate the non-relativistic mean-field potential U™' in the local rest frame 
(LRF) of nuclear matter which is defined by the frame of reference with vanishing local 
vector baryon current {j(r,t) = 0). Discarding vector potentials in the LRF we then 
equate the expressions for the single-particle energies using the non-relativistic potential 
and the scalar potential U by 

^Jp^+m^ + [/"'■(f, p) = \Ip^ + {m + U {r,p)f. (11) 

Eq. (|Tl] ) now allows to extract the scalar mean-field potential U (r, p) which we will use 
throughout our calculations for the baryons. Due to the relativistic dispersion relation 



for the quasiparticle (11) the scalar potential U{f, p) now has definite Lorentz-properties 
contrary to This enables us to guarantee energy conservation in each two-body 

collision and in resonance decays (eq. Ni + N2 ^ + N4) as 

El + E2 = E3 + £^4 



pf + ml^ + y p^ + m^^ = \lpi + rnY + ^Jpl + mY. (12) 

Since the collision integrals (P) are evaluated in the center-of-mass system of the colliding 
particles or in the restframe of the decaying resonance, while the testparticle equations 
of motion are generally integrated in the center-of-mass system of the heavy- ion reaction, 
energy conservation is not a priori fulfilled when employing non-relativistic potentials 
without definite properties under Lorentz transformations. 

For our calculations we will use a (momentum-dependent) equation of state (EOS) for 
nuclear matter with an incompressibility of K = 250 MeV ( A = -29.253 MeV, B = 57.248 
MeV, C = -63.516 MeV, r = 1.760, A = 2.13 1/fm.) The actual calculation of the mean- 
field potential according to eq. ([lO|), however, is too involved for practical purposes since 
one has to perform double integrations for the momentum-dependent part of ?7"^(r, p). We 
thus determine the potentials within the Local Thomas-Fermi approximation, where the 
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integral over the phase-space distribution function in (10) can be performed analytically 



[2|. With 



where the factor 4 stems from the sumation over spin and isospin, the integral for the 
momentum-dependent part of the potential f/"*" reads: 



.3^/ fi^^P') 4 3 
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1 + 



: i''n 7^ lo 

2pK (p_p^)2+A2 ^ ' 



2pF 

A 
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We employ a smeared baryon density |^] when evaluating expression (13) rather than that 
obtained directly from the test-particle distribution in order to avoid unphysical statistical 
fluctuations in the density and in the mean-field potentials. 

3.2 The Coulomb Potential 

Charged baryons and mesons are additionally exposed to the Coulomb potential 14 (r) 
generated by all charged particles. Since in our present approach mesons are propagated 
as free particles with respect to the nuclear mean-field, the Coulomb force 

F,{r) = -qVrVc{f^ (14) 

is the only force acting on a meson with charge q. For charged baryons the force Fc ([l^ 
represents an additional term in the equations of motion (^). The Coulomb potential 
Vc (r) is obtained by solving the Poisson-equation by means of the Alternating- Direction 



Implicit Iterative ( AD I-) algorithm |23]. 



3.3 Resonance Properties and Decay Widths 

Within the CBUU-approach the resonances are treated as "on-shell" particles with respect 
to their propagation and the evaluation of cross sections. To account for their "off-shell" 
behaviour, we distribute the resonance masses according to a Lorentzian distribution func- 
tion (see sec. (p.5|)), which is determined by the mean resonance masses Mr and the total 
and partial decay widths Tn at mass M/j. For the baryonic resonances these explicit pa- 
rameters are given in table ||. For the p-meson we use a mean mass Mr = 770 MeV/c^ 
and for the decay width at resonance Tpi = 118 MeV. The corresponding values for the 
cT-meson are Mr = 800 MeV/c^ and Tr = 800 MeV. 

In the following we list the parametrizations used for the decay widths T{M). 

• Ivr-decay width for the A(1232) 

For the A(1232)-decay we adopt the parametrization given by Koch et al. p3] 



Moniz\ 




i<i) = ^R^[-] r-5r^\ > (15) 



where M is the actual mass of the A (1232) and Ma = 1232 MeV/c^. q and qr are 
the pion three-momenta in the restframe of the resonance with mass M and Ma, 
respectively. The parameter 5 in the cutoff function has a value 5 = 0.3 GeV/c^. 
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Table 1: Decay widths and decay channels for the baryonic resonances. The data for 
these channels are taken from and [^] in case of the A^(1535). The second column 
contains the averaged matrix-elements for the production of the baryon resonances in 
nucleon-nucleon collisions (cf. sec. 3.5.2| ). 



Ivr-decay width for the higher baryon resonances 

The l7r/7/-decay widths for the higher baryon resonances are given by 



where I is the angular momentum of the emitted pion or rj and q and qr are the 
momenta of the pion or in the restframe of the decaying resonance as defined 
above. In this case, we use 

5^ = {Mr -Mn- m^f + ^ (17) 

27r-decay width for baryon resonances 

The 27r-decay of the higher baryon resonances is described by a two-step process. 
First, a higher lying baryonic resonance decays into a A(1232) or A^(1440) and a pion 
or into a nucleon and p- or cr-meson. The new resonances then propagate through 
the nuclear medium and eventually decay into a nucleon and pion or a nucleon and 
two pions 

rb ^ Ntttt. (18) 

Here R denotes the higher baryonic resonances, r stands for a A(1232), A^(1440), p or 
a and b is a nucleon or pion, respectively. Due to the fact that a further resonance 
appears in the first step of reaction ( p^sj) one has to integrate the corresponding 
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distribution function over the mass fx of the intermediate resonance r in order to 
obtain the 27r-decay width 

(19) 



'{MR-MN-2m^f + 5^V 



\{M-Mn- 2m^y + 52 

where is the branching ratio for the 27r-decay of the baryonic resonance R and 
Pf denotes the momentum of r and b in the restframe of R. Since the integral in 
eq. (^) for high resonance masses M is proportional to M, we introduce a cutoff 
function to avoid T^^raiM) to diverge for high masses. The parameter 6 used is 
0.3 GeV. 

decay width for meson resonances 

The decay width of the meson resonances is parametrized similarly to that of the 
A(1232), 

Mr f q\^^^+^ q^r+S^ 



where and M are the mean mass and the actual mass of the meson resonance. 
q and qr are defined as in eq. ( [T5| ) while Jr is the spin of the resonance and T^ 
the decay width for a resonance with mass Mr- For the parameter 6 in the cutoff 
function we use again 6 = 0.3 GeV. 

3.4 Meson-Baryon Cross Sections 

In order to describe meson-baryon scattering in the framework of our resonance picture 
we use a Breit-Wigner formulation for the cross sections (eq. vrA^ irN) 

_ 2Jr + 1 47r S TR^ab ^R^cd 

aab^R^cd - (25^ + i)(2S, + 1) p2 (, _ M2 )2 + s TL • ^ ^ 

In eq. (|2l|) ab and cd denote the baryon and the meson in the final and initial state of the 
reaction and R is the intermediate baryon resonance. Jr, Sa and Sb are the spins of the 
baryon resonance and the particles in the initial state of the reaction. 
The solid line in fig. |^ shows the total vr" — p-cross section within our model in comparison 
to the experimental data from p^. To calculate this cross section we replace the partial 



widths T R_^cd in eq. (21) by the total widths of the baryonic resonances and sum up the 
contributions from all resonances. The dashed, the dotted and the dash-dotted lines in 
fig. |l| show the contributions from the A(1232), the A^(1440) and the A^(1535) separately. 
In fig. ^ we display the resulting cross section for the Ti~p rjn reaction. Here only the 
A^(1535) contributes since this is the only baryonic resonance that couples to the rj within 
our model space (cf. table |l|). Both cross sections are obviously well described up to k, 
1.0 GeV. 

Eq. (^) is also used to determine the cross sections for a- and p-production in vr — vr- 
collisions weighting with the corresponding spins of the mesonic resonances and the pions 
in the initial state of the reaction. 
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Figure 1: The calculated total vr — p-cross section in comparison to the data from |26] 
is given by the solid line. The dashed, the dotted and the dashed-dotted lines show the 
contributions from the A(1232), A^(1440) and iV(1535), respectively. 
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3.5 Inelastic Baryon-Baryon Cross Sections 

In this section we describe the concepts and parametrizations used to implement the cross 
sections for resonance, pion and rj production in the CBUU-model. 

3.5.1 The NN iVA(1232) Cross Section 

For the NN A^A(1232) cross section we employ the result of the OBE-model calculation 
by Dimitriev and Sushkov ||l^ using u- and t-channel Born-diagrams. The parameters of 
the model (NNn-, A^A(1232)7r-coupling constant and a parameter in the vr-formfactor) are 
chosen to reproduce the experimental pp A^A"*"^ cross section. For the implementation 
of this cross section into the CBUU-model we replace the parametrization for the A(1232)- 
width given in |jl^ by the Moniz parametrization ([l5|). An example for the resulting mass 
and angular differential cross sections for an invariant energy of 2.31 GeV is given in fig. 
|3| in comparison to the experimental data ||l^. The cross sections for the other isospin 
channels follow from the pp nA++ cross section by applying isospin symmetry. 

3.5.2 The NN NR Cross Section 

In order to obtain the production cross sections for the higher baryon resonances in 
nucelon-nucelon collisions we fit the corresponding matrix-elements to available data for 
r]-, Ivr-, p- and 27r-production in nucleon-nucleon reactions. Therefore we assume that 
the 7]-, Ivr-, p- and 27r-production in nucleon-nucleon collisions above the A(1232) exci- 
tation proceed only through these resonances via subsequent two-step processes. The r], 
single pion and p production are described by the creation of a baryonic resonance in a 
nucleon-nucelon collision and its subsequent decay into r], pion or p 

NN ^ NR^ NN-n/'q/p. (22) 

We assume that the 27r-production in nucleon-nucleon collisions proceeds either through 
the excitation of two A(1232) and their subsequent decay into a nucleon and a pion or 
through the excitation of higher lying baryonic resonances and their subsequent decay into 
a nucleon and two pions (see section ( |3.3D ). 



The general expression for the cross section 

For the derivation of the cross section we assume that in a collision of two nucleons 
a+h — > R+c a baryonic resonance R and a nucleon c are produced. Then resonance R 
decays into a two-body final statei? de. Assuming spinless particles the invariant 
matrix element is given by 

M — -Mab-^Re Pr M-R^cd, (23) 

where Pr is the propagator of the intermediate baryonic resonance R and Mab-^Re 
and Mpi^cd ^-re the matrix elements for the reactions a + 6 — > i? + e and i? — > c + d, 
respectively. We start from the general expression for the cross section 

dcrab-,cde = T— ^ ^^{Pa + Pb - Pc - Pd - Pe)\Mab^cde\'^ 
4: Pi VS 



d^Pc d^pd d^pe 
(27r)3 2Ec (27r)3 2Ed {2-K f 2Ee 



(24) 
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Figure 3: Comparison of the pp — > nA(1232) angular and mass differential cross section 
at an invariant energy of 2.31 GeV (solid lines) with the experimental data |l^] (filled 
circles) . 
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where -y/i is the invariant energy of the particles in the initial state and pi is their 
CMS momentum. \M.ab^cde\^ is the square of the invariant matrix element averaged 
over the spin of particles a and h and summed over the spins of the particles in the 
final state of the reaction. Assuming that the square of the matrix element factorizes 
as 

\Mab-,cde? = \Mab^Re? \Pr? \MR^cd\\ (25) 

we obtain for the cross section as a function of the resonance mass /x 

= '^--^^(^^ n ^,^-MlY + ,^Ti; 

Here, (Tab^Re{^^) is the cross section for producing a baryonic resonance R with a 
fixed mass /i 

(rah^Reifj) = -^tA / dVLpf\Mab->Re{lJ')?, (27) 

64 TT^ S ^Ji J 

with pf and pi denoting the CMS momenta of the particles in the final and initial 
state of the reaction a + 6 — > i? + e, respectively, while s is the squared invariant 
energy of this reaction. When evaluating the production cross sections for the higher 
baryonic resonances we assume \M.ab->Re{lA\^ to be constant for all baryonic reso- 
nances except for the A(1232) (cf. section (|3.5.lD ). 



//-production cross section 

Since in our model only the A^(1535) couples to the r^-meson, we obtain the cross 
section for r/-production in nucleon-nucleon collisions from eq. ( p6|) by using the 
iV(1535) as the intermediate baryonic resonance and regarding NN-q as the final 

state cde. The unknown squares of the matrix elements \M.ab~>Re\^ are obtained by 
fitting the available experimental data for ry-production in nucleon-nucleon collisions. 
The matrix element for A^(1535) production in proton-proton collisions then gives 



M 



pp^pAf+(1535) 



167r X SmbCeV^ 



(28) 



The result of this fit for the reaction pp pprj is displayed in fig. ^ in comparison 
to the data |2?, It is known from |27] that the cross section for ?7-production 
in proton-neutron collisions is about a factor 5 larger than that for proton-proton 
collisions. Consequently we use 



M 



pn-»pAf(1535) 



M 



pn^n.Af+(1535) 



167r X 40mbGeV^ 



(29) 



Ivr-production cross section 

In order to evaluate the Ivr-production cross sections in nucleon-nucleon collisions 
we sum up the contributions from all resonances contributing to a specific channel 
incoherently. For A(1232)-production we use the cross section given in sec. ( ^.5.1 ) 



and for the other baryonic resonances we use eq. (26) integrated over the resonance 
mass with a and b denoting the two nucleons in the initial state and c, d and e 
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Figure 4: Cross section for the reaction pp pprj obtained from the resonance model 
using the matrix element given in eq. (p3) (solid line) in comparison to the experimental 



data 1^ 



the NNtt final state. Introducing the proper isospin coefficients we obtain 



with 

0"3/2 





-^p p tt" 


2 1 

30-3/2 + 3^1/2 


^pp~ 


+pn7r+ 


10 2 

^0-3/2 + 3^1/2 


^p n— 


^ppTT~ 


1 1 

3^3/2 + 3^1/2 


^p n- 


-^pn-K^ 


4 1 

3^^3/2 + 3^1/2' 


R+-- 




-(1232) 



(30) 
(31) 
(32) 
(33) 



1 \M\' r ■ , V- 

H ■ ■ — / "MPf X / o 

4 IQirpi s JM^+m^ ^^f^,3,) (^2 _ M^f + ^2 ^J^) 



(34) 



and 



0-1/2 = PP^PR^ 
Ir=\ 



Wttp.s JM^+m^ ^ V {^J?-Mlf + ^J? r| (/x) 



(35) 



, _ 1 



There is no explicit factor r^y-Tr in the contribution of the A(1232) to 173/2 (|34| ) be- 
cause it decays with 100% probability into a nucleon and a pion. 
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Figure 5: Cross section for the reaction pp ppp^ obtained from the resonance model 
using the matrix elements given in table |^ (solid line) in comparison to the experimental 
data p|. 



p-production cross section 

The invariant matrix-elements for the production of baryon resonances which can 
decay into a nucleon and a p (cf. table ^) are obtained by a fit to the experimental 
data for the reaction 

pp — > ppp^. (36) 

Similar to eq. (30) we write the p'^-production cross section as a sum of contributions 
of I = 3/2- and / = 1/2-resonances, 

2 1 

where CT3/2 and are defined as in eqs. ( ^4|) and (|35|) , while the sum extends only 
over resonances with a non-vanishing branching ratio for the decay into a nucleon 
and a p. The resulting matrix-elements are quoted in the second column of table |^. 
The corresponding p'^-production cross section is shown in fig. S in comparison to 



the experimental data |26]. 



The 27r-production cross section 



As already stated in sec. (3^) the 27r-production cross section in nucleon-nucleon 
collisions is described via the excitation of higher lying baryonic resonances and 
their subsequent decay into a nucleon and two pions according to the branching 
ratios given in table |l| and the corresponding decay widths: 

1. NN ^ NR^ iVA(1232)7r ^ NNtttt (38) 

2. NN ^ NR^ iViV(1440)7r ^ NNtttt (39) 

3. NN ^ NR^ NNp NNtttt (40) 

4. NN ^ NR^ NNa NNtttt. (41) 
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PP ^ PP" " 




TTn if TTn Tr"^ 


Till ^ T>7T7r^~'7r 

jJIL ' jJILn i\ 


PP ^ PPT^ ^ 




o/y 


z/y 


2/9 


o/y 


z/y 


Hn 


4/q 


1/y 


4/9 




4/9 


ns 


1/3 





2/3 


1/3 


2/3 




2/3 


1/3 





2/3 





di 


26/45 


2/45 


22/9 


52/45 


17/45 


d2 


2/9 


2/9 


14/9 


4/9 


5/9 


d3 


2/3 





10/3 


2/3 


1/3 


di 


















Table 2: The products of the isospin coeficients for the three steps of the reactions ( ^8|) - 
( |4l| ) for those isospin channels where experimental data are available. 

Here R stands for the higher lying baryon resonances. In addition to these processes 
we take into account the 27r-production via the excitation of two A(1232), 

NN A(1232)A(1232) ^ A^iVvrvr, (42) 



adopting the cross sections from ref. [1^. For the description of the 27r-production 
cross section in nucleon-nucleon collisions we define cr^^ ^1/2 ~ sim- 



ilarly to eqs. ( p^ and (35) replacing Tr^attt by the corresponding decay width Tj 
{i = 1, ..,4) responsible for the second step of the reactions (|38|) to (p|). The cross 
sections then read 

4 4 

(^NN^NNtt-k = ^ni<^l/2 + ^dia\i2 + (^NN-^ A(1232)Ail232)^NNmr, (43) 
i=l i=l 

where the rii and di are the products of the isospin coefficients for the three steps of 
the reactions (^) - (pi]). The corresponding factors are listed in table |2|. 

Using the matrix elements already determined by the fits to the t]- and /9-production data 
we now adjust the unknown matrix elements IMnn^nrI"^ to reproduce the cross sections 
for Ivr- and 27r-production in nucleon-nucleon collisions. The resulting Ivr-cross sections 
are shown in fig. ^ (solid line), where the contributions from the A(1232) (dashed line), 
the sum of all contributions from the isospin-1/2 resonances (dotted line) and the sum 
of all higher isospin-3/2 resonances (dash-dotted line) are displayed separately. Evidently 
the pion-cross sections are fitted in our resonance model reasonably well up to invariant 
energies of 5 GeV. 

In fig. 1^ we plot the resulting 27r-production cross section ( ^3| ) for the isospin channels 
where experimental data are available |2^. As in case of the Ivr-production channels 
the 27r-data can be reproduced well within our multi-resonance model using the matrix 
elements given in table m. 

3.5.3 The NR NN Cross Section 

The cross sections for the reaction NR — > A^A^, where R stands for the A(1232) as well as a 
higher baryonic resonances, are given by an expression corresponding to eq. (^) for a two- 
body final state. For the matrix elements involved we adopt those from the OBE-model of 
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Figure 6: The fitted Ivr-production cross sections (solid line) for different isospin channels 
in comparison to the data |26|; contributions from the A(1232) (dashed), the isospin 1/2 
resonances (dotted), the higher isospin 3/2 resonances (dashed-dotted) . 



Dimitriev |]T5| for the A(1232) and those from the simple resonance model (sec. (3.5.2)) 
in case of higher baryon resonances. Thus, in comparison to previous implementations 
of the BUU-model 



25, 



we do not have to employ a detailed balance prescription to 
deduce the NR — > A^A^-cross section, since the matrix elements are known as a function 
of the invariant energy and the resonance mass. 



3.5.4 The NR NR'-Cross Section 

For the collisions of a nucleon and a resonance R leading to a nucleon and a different 
resonance R' we use for all resonances (including the A(1232)) the average of the matrix 
elements obtained for the reactions A^A^ NR and A^A'^ —>■ NR' . In analogy to eq. (|2^ ) 
this leads to the following expression 

^ 0.5( Aljviv^iVii +\M 

NN^NR' 

)2(2J^,+1) 

aNR^NR' = I 167^-^ (44) 

where / accounts for the proper isospin coefficients (cf. table |3|) and stands for the 
spin of the resonance in the final channel. 



3.5.5 The A^A^ NNn-Cross Section 

The A^A^ A^i?-cross sections are determined by fitting the experimental data for Itt- 
production in nucleon-nucleon collisions above the A(1232). As compared to the data the 
resulting Ivr-cross sections are slightly too low just above the vr-production threshold (cf. 
fig. |8|). In order to compensate for this deficit we attribute the difference between the 
data and the cross sections obtained from the resonance model to the cross section for 
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Figure 7: Cross sections for 27r-production in nucleon-nucleon collisions for different isospin 
channels obtained within the simple resonance model in comparison to the experimental 
data @. 
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A^ 


1/4 


N+ A++ 
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1 


iV+ A+ 


AftJ A++ 


3/8 


iV+ A+ 


N+ A+ 


5/8 


N+ A'J 


A^ 


1/2 


iV+ A- 


N+ A- 


5/8 



Table 3: Isospin coefficients / for baryon-baryon collisions. N and A represent isospin 
1/2 and isospin 3/2 particles, respectively. The isospin channels not listed explicitely are 
determined by assuming isospin symmetry. 



channel 


A [mb] 


a 


b 


ni 


"2 


p p ^ p p vr*^ 


61.3 


1.52 


2.50 


6.18 


3.48 


p p ^ p n TT^ 


122.6 


1.52 


2.50 


6.18 


3.48 


p n ^ p p TT^ 


24.9 


3.30 


0.85 


1.93 


0.002 


p n ^ p n TT^ 


7.25 


0.88 





2.31 


3.64 



Table 4: Parameters obtained for the direct vr-production cross section n^n N n (eq. 



direct (s-wave) vr-production in nucleon-nucleon collisions (A^A^ NNn). The resulting 
difference (s-wave) cross section can be fitted by the expression 

aNN^NNA^) = Ax^' e-^-^"'+'^\ (45) 

with 

^/s — 2Mm — rUn 
5G^V ' 

using the parameters given in table ^. Now adding ajy n incoherently to the Ivr- 
production cross section from the baryon resonance decays then yields the total cross 
section depicted in fig. |8| by the solid line which now also gives a fit at threshold. 



3.5.6 The A^A^vr A^A^-Rate 

Assuming that the invariant matrix element for the reaction A^i + A^2 ^3 + A4 -|- vr 
depends only on the invariant energy one can write the cross section (JNiN2^N3N4,-k (see 
sec. ( 3.5.5| )) in the following form |13] 

64(2 7r)3pi^'' J 

Here S]y-^^^]y^ is the symmetry factor for A^i and A'2, Pi is the CMS momentum of the 
nucleons in the initial state and 
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Figure 8: The Ivr-production cross section taking into account only the contributions due 
to the decay of baryonic resonances (dashed hne) and the total cross sections obtained by 
adding the direct pion production in nucleon-nucleon collisions (solid line) in comparison 
to the experimental data |2(:]. 



The transition rate for a pion being absorbed by N'^, and A'4 is given by 

^4 (mi +pjv,2 -m3 -m4 -ptt) \m.\^ 

Wfi = {2tt) , (47) 



where the normalization volume V contains 2E particles |30|. Multiplying this equation 
with the phase-space factors for the two nucleons (iVi, A'^2) in the final state and taking 
into account that the pion reacts with two nucleons from the surrounding nuclear medium 
gives 

V^N,N,^N,N, - Sm,,n,j^\M\ (48) 
where pM2, ^-nd pN i are the corresponding local neutron or proton densities. 



3.6 Elastic Baryon-Baryon Cross Sections 

For the elastic nucleon-nucleon cross section we use the conventional Cugnon parametriza- 
tion [|, m 

OMN^NN = ^ ^ X,,,^ + 20 1 mb. (49) 

\ ^ + GcV 

The cross sections for nucleon-resonance scattering (with the same baryonic resonance in 
the initial and final channel) are determined in the following way: assuming an isotropic 
angular dependence for ( ^9|) we fix the matrix element squared by 



\Mnn^nn\^ = IQt^s\ ^— — + 20|mb. (50) 



35 



1 + 



v/i-2Af 
GeV 



JV 



20 



Now inserting (|50|) in eq. ( ^61) we obtain for the elastic nucleon-baryon scattering cross 
section 



\Mnn^nn\ f ^ 2 f^^Tni^) 



fd^p/^ ™ . (51) 



We note that the cross section (^) in addition to elastic scattering also allows for a change 
in mass of the resonance in the scattering process. 

4 Results for Nucleus-Nucleus Collisions and Comparison 
to Data 

As described above the current implementation of the CBUU model takes into account 
baryonic resonances up to masses of 1.95 GeV/c^ for the first time. Hence, we start 
our discussion with an analysis of the effects of these resonances on the resulting pion 
spectra. Since the higher baryon resonances decay ultimately into a nucleon and one 
pion or into a nucleon and two pions (cf. sec. ( |3.5.2|) ) we expect the contributions from 
those resonances to show up in the low energy as well as in the high energy regime of 
the pion spectra. An enhancement of the pion yield in the low energy regime is expected 
due to the 27r-channels {NN A(1232)A(1232) NNtttt, NN NR ^ Ntttt, sec. 
( 3.5.2[) ). Pions produced via the decay of two A(1232), that are excited simultaneously 



in a nucleon- nucleon collisions, or via the 27r-decay of the higher resonances have a lower 
momentum than those stemming from the reaction A^A^ — > NR NNtt. Thus we expect 
an enhancement of the low energy pion yield. On the other hand, once a higher baryon 
resonance is excited and subsequently decays into a nucleon and a single pion, this pion 
will have a higher momentum in the restframe of the resonance than those emitted in 
A(1232)-decays. 

In fig. ^ we show the effect of the 27r-channels on the low energy regime of the 
pion spectra. Fig. |9| displays the total vr- multiplicity weighted with 1/pt as a function 
of the pion transverse momentum px for a central Au + Au reaction at 1.5 AGeV. The 
solid line represents the result of a full CBUU-calculation and the dashed line is obtained 
by switching off the 27r-channels. This is done by neglecting the cross section AA^ — > 
A(1232)A(1232) as well as that for the time reversed reaction. The effect of the 27r-decay 
of the higher resonances is eliminated by allowing only the decay into a nucleon and a single 
pion, i. e. (TiT^/Ttot = 1 and T2tt = 0). In addition the matrix elements for resonance 
production in nucleon-nucleon scattering have been refitted in order to guarantee the 
correct description of the Ivr-production cross sections in nucleon-nucleon collisions. The 
resulting value for the matrix element squared in the latter case is 



\MN,N^N,Rf = 167r X 9mbGeV^. 

As can be seen from fig. ^ at 1.5 GeV/A the 27r-channels increase the pion yield at low pT 
by 20 — 25%. The effect vanishes already at a transverse momentum of 0.25 to 0.3 GeV/c. 
For beam energies about 1.0 GeV/A we find that this effect is reduced to ~ 10%. After 
investigating the effect of the 27r-channels on the vr-spectra we now discuss the impact of 
the higher baryonic resonances at high pion kinetic energies. In this respect we display in 
fig. |l^ the 7r~ -multiplicity in central La + La collisions at 1.35 AGeV as a function of the 
center-of-mass kinetic energy of the pions at 80° < QcAis ^ 120° obtained from a CBUU- 



calculation in comparison to the experimental data from |31|. We are able to reproduce 
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Figure 9: Total pion multiplicity weighted with 1/pt as a function of the transverse 
momentum pT for a central (b = 0.0 fm) Au + Au collision at 1.5 GeV/A. The solid line 
shows the result of the full CBUU-calculation and the dashed line corresponds to the result 
when omitting the 27r-production channels as described in the text. 

the experimental spectrum; especially the slope at high pion energies {Tcms > 0.5 GeV) is 
described well by the calculation. That this slope is influenced significantly by the higher 
lying resonances is shown in fig. |ll|, where we display the calculated pion multiplicity for 
a central La + La collision at 1.35 GeV/ A as a function of the center-of-mass pion kinetic 
energy. The upper solid line corresponds to the full calculation while the dotted line dis- 
plays final pions stemming from direct pion production in nucleon-nucleon collisions or 
from A(1232)-, p- or u-decay. The dashed line corresponds to the sum of the final pions 
stemming from the decay of higher baryonic resonances; the histograms show individually 
the contributions from single resonances to the pion yield: A^(1440 (solid line), A^(1535) 
dashed line, A/^(1520) dash-dotted line and A(1600) dashed line. The sum of the contri- 
butions from higher resonances for low kinetic energies is in the order of 10 — 20% while 
the yield above Tcms ~ 0.6 GeV is fully dominated by pions originating from the decay 
of higher resonances. 

We now turn to a comparison of our calculations to the data on 7r~-production in 
Ar + KCL collisions at 1.8 GeV/A obtained at the BEVALAC. All data shown in the 
following are taken from ref. [|l^]. In the experimental analysis a central and a minimum 
bias event class have been used. In order to be able to relate the CBUU-results to the 
experimental events we determine these event classes by comparing to the -xp-spectrum. 
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Figure 10: The vr" -multiplicity for central La + La collisions at 1.35 GeV/A as a function 
of the center of mass pion kinetic energy for 80° < @CMS ^ 120°. The solid histogram 
displays the result of the calculation while the shaded areas indicate the statistical error 
of the calculation. The data (squares) are taken from [31|. 



The inclusive cross section is obtained by the integration of the multiplicities n(b) over 
the impact parameter b via 

da _ dni}^^^^^ ^^^^ 



dpT Jo dpT 

In order to identify the "central event" class we thus integrate up to a maximal impact 



parameter bmaxi which we determine by fitting the data. Fig. 12 shows these cross 
sections as a function of the transverse pion momentum pT from the CBUU-model for 
bmax from 1.6 to 2.3 fm. In this ^max-regime one can see that the shape of the calculated 
spectrum does not depend on the events used; one only observes a shift of the spectrum 
in magnitude when including more events. Since the integrated spectrum is reproduced 
well for bmax = 2.1 fm, we will use events with b < 2.1 fm for the comparison of the 
CBUU-calculations with the data. 



From fig. 12 we see that the spectral form of the data is described well over the 



whole momentum range except for the transverse momenta between 0.3 and 0.55 GeV/c 
where the data are overestimated. This overall agreement is also found when looking at 
as a function of the pion kinetic energy in the CMS (cf. fig. |l^), although we find 
that the distribution resulting from the calculation seems to be shifted by ~ 25 MeV to 
higher energies. In fig. 14, furthermore, we show a double differential spectrum where the 
data are indicated by the squares and the solid histrogram corresponds to the result of 
the CBUU-calculation obtained by scaling the original calculation in units mb/GeV/sr by 
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Figure 11: The 7r~-multiplicity for central La + La collisions at 1.35 AGeV as a function of 
the center-of-mass pion kinetic energy. Solid line: CBUU-calculation, dotted line: fraction 
of pions stemming from A(1232)-, p-, fi-decay or from direct pion production in nucleon- 
nucleon collisions; dashed line: fraction of pions originating from decays of higher baryonic 
resonances. Histograms: fraction of pions from specific resonance decays: solid A^(1440), 
dotted 7V(1535), dash-dotted 7V(1520) and dashed A(1600). 



a factor 1.9. Again we observe a good agreement for the form of the spectrum, although 
again the calculation seems to overestimate the data between 0.3 and 0.5 GeV. 

After convincing ourselves that the CBUU-model reproduces the spectral form for the 
7r~ -distributions rather well we now turn to differential pion angular distributions. In fig. 



15| we display J^^n for the inclusive vr -yield, i.e. integrating eq. (^) up to oo (top), and 



for the central event class (bottom) in comparison to the data ||l^. The circles represent 
the data while the squares indicate the results of the calculation. The solid lines are fits 
to the calculations employing the functional form [^] 

const X (1.0 ^acos^e). (53) 



dcos9 



For both event classes the data are reproduced remarkably well. When considering the 
minimum bias event class we find for the anisotropy parameter a ~ 1.03, while a is reduced 
to ~ 0.51 when looking only at central events. As stated in [|^] this decrease of anisotropy 
- when going from minimum bias to central events - can be understood as an effect of the 
centrality of the heavy-ion collision because for minimum bias events semi-peripheral and 
peripheral collisions are weighted stronger than central collisions (c.f. eq. ([5^)). Pions 
produced in semi-peripheral and peripheral collsions are more likely to originate from first 
chance A^A^ — > AA NNtt collisions than those produced in central collisions. Thus 
the anisotropy introduced in these first chance collisions prevails over the more isotropic 
pion distributions originating from central collisions. This line of argument is supported 
by the fact that the anisotropy coefficient is reduced to half of its value when considering 
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Figure 12: The 7r~ -cross section for central event classes as a function of the transverse 
momentum pT- The squares indicate the data from |1C] while the lines show the result of 
a CBUU-calculation with different bmax '■ solid line bmax = 2.1 fm , dashed line bmax = 1-7 
fm, dotted line bmax = 1-9 fm, dash-dotted line bmax = 2.3 fm. 
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Figure 13: The differential cross section E-^^ for the central event class {bmax = 2.1 fm) 
as a function of the pion-kinetic energy. The squares indicate the data from |l^ and the 
solid line corresponds to the calculation. 
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Figure 14: The double differential vr'-cross section for the central event class as a function 
of the pion-energy. The squares represent the data, while the solid line corresponds to 
the result of the CBUU-calculation. The shaded areas indicate the statistical errors of the 
calculation. 



only central collisions (see. bottom of fig. 15). In central heavy-ion collisions a high 
density regime is formed where pions and baryonic resonances are produced and absorbed 
repeatedly. Thus, pions from central collisions result mainly from multi-step processes 
and hence are emitted more isotropically than pions from peripheral collisions. 
Finally, we look (for the central event class) at the energy dependence of the anisotropy 

m 



parameter a assuming as m 



a{E) ~ 1 + a{E) cos^O. 



(54) 



In fig. 16 we show a{E) as a function of the pion-kinetic energy. The result obained 



from the CBUU-model is depicted by the squares while the circles indicate the data. 
The overall functional form of the data is reproduced by the calculation except for the 
pronounced peak for pion kinetic energies of 0.25 to 0.35 GeV. The CBUU-model also 
shows the increase in anisotropy for energies up to 0.3 GeV from zero to an a{E) of ~ 1.0. 
Above these energies we observe a slight decrease of the anisotropy coefficient a(E). 
In fig. ^ we finally compare our calculations to the more recent data on vr'^-production in 
Au + Au collisions at 1.0 AGeV from the TAPS-collaboration. The solid histogram shows 
the result of our calculation, while the shaded areas indicate the statistical error of the 
calculation. The squares represent the data form |32, 5^. Here we find that our calculation 



is in good agreement with the experimental data, only slightly overestimating the spectrum 
in the region of pj- ~ 0.2 GeV/c. In comparison to earlier calculations within the BUU- 
model of ref. (sj] the low p^-behaviour of the pion-spectrum has improved. In addition 
we are now able to reproduce the high p^^-data while the earlier calculation overestimated 
the high p^-spectrum. The results obtained by Bass et al. within the IQMD-model [0] 
are ~ 20 % higher than ours in the low p^-region, while both calculations are well in line 
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Figure 15: The number of pions as a function of cosO in the CMS for minimum bias (top) 
and central (bottom) events. The data are represented by the circles and the results of 
the CBUU-calculations are given by the squares. The solid lines represent fits of the form 
(l5l) to the CBUU-calculations. 
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Figure 16: The anisotropy parameter a (see text) as a function of the pion-kinetic energy 
for the central event class. The circles represent the data and the squares indicate the 
result obtained by fitting the CBUU-results in the regime 0.0 < cos9 < 0.9 by const{l + 
a{E)cos^6) for fixed pion kinetic energy. 

for the other p^^-regions. 

5 Summary 

In this paper we have presented a new implementation of the CBUU-transport-model 
for the description of relativistic heavy-ion collisions up to energies of about 2 GeV/A 
including baryonic resonances up to masses of 1.95 GeV/c^ as weU as 7r,r7 and p mesons. 
In addition to earlier implementations we also include 27r production channels via NN 
AA — NNtttt and the 27r-decays of higher baryon resonances. The inclusive Itt and 27r 
cross sections from nucleon-nucleon collisions are found to be well reproduced within our 
multi-resonance approach. 

A detailed analysis of nucleus-nucleus collisions from 1-2 GeV/A shows that the 
27r production channels increase the pion yield essentially in the region of low transverse 
momenta pT while the one pion decays of the higher baryon resonances are dominating 
the high pT or energetic part of the pion spectrum. Thus the energetic part of the pion 
spectrum directly reflects the abundancy of higher baryon resonances in the compressed 
stage of the nucleus-nucleus collision. 

In comparison to the Ar + KCl-data at 1.8 AGcV for vr" -production obtained at 
the BEVALAC we found an overall agreement of the spectral 7r~-distributions from the 
CBUU-model with the data except for the low p^-regime where the data are slightly 
underestimated. In comparing our calculations with the 7r°-data of the TAPS-collaboration 
for Au + Au at 1 GeV/A we find that the CBUU-model is in good agreement with the 
experimental spectrum, except for transverse momenta of 0.2 GeV/c, where the data 
are overestimated. Furthermore, we have shown that the experimentally observed angular 
anisotropics for pions are reproduced well and can be understood in the framework of the 
multi-resonance picture. 
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Figure 17: Cross section for vr'^-production in Au + Au collisions at 1.0 GeV/A for 0.52 < 
Yiab ^ 0.84 as a function of the pion transverse momentum px- The solid histogram 
displays the result of the CBUU-calculation, while the shaded areas indicate the statistical 
error of the calculation. The data TAPS-data are represented by the squares. These data 
from ref. |^2[ have been renormalized by a factor 0.6 determined by a recent new analysis 
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The NN NR-cToss section 



reaction: (e.g. NN NN{1U0) NNn) 

a -\- b ^ R -\- e c + d + e 
R ^ c + d 



matrix element for spinless particles: 

M M PrM 

Pr. propagator of intermediate resonance 



general expression for the cross section 



dUah^cde = -^—^S'^{Pa+Pb-Pc-Pd-Pe)\Mab^cdef 

d^Pc d^pd d^pe 

X ■ 



[2tiY2EA2tiY2EA2tiY2E, 



|A^a6_^crfe| '■ square of the invariant matrix element averaged over the 
spin of particles a and h and summed over the spins of the particles 
in the initial and final state of the reaction. 



1 



assume: factorization of square of matrix element: 



\M \ \Pr\ \M 



mass differential cross section 



do'ah-^Re-^cde / \ 2 /i^ T 



/ \ - ^ ^ R^cd 

<^ab^Re[l^) 



2 

tot 



with the cross section for resonance production 



1 I- 2 

aah^Reip) = -TT—. / dVt p f\Mah->Re{p)\ 

64 TT^ spi J 



Pi, pf initial and final CMS momenta, ^/s invariant energy of the re- 
action 



Assume \Mab->Re{l^)f = const. 



fit to experimental cross section for Itt-, 27r-, rj- and /)-production 



The NN NNrj-cToss section 



fitting matrix element to pp pprj data yields: 



M 



= IGtt X SmbGeV^ 



experimentally: //-production in pn-collisons 5 times larger than in 
j9j9-collisions 



IGtt X 40mbGeV2 



The NN NNtv-ctoss section 



Incoherent sum over the contributions of all resonances 

define: 



T —3 



and 



1 \Mf rVs-MN 
0"j9j9^j9A+(1232) + TT7 L ^ 

iJ^A(1232) (/i2 - M|) + /i2 rfj ^,^(/i) 



(^1/2 = VP ^pR'^ 

T _1 



assume incoherent summation of contributing isospin channels 
e.g.: 

[pn\Hi,t\N*^n)\^ = 

' 1 



^ < 101 + 4^ < 00|Wi„t (^110 > 



,72 
4 + 4 



72 



,72 



72 



1 , ^2 

00 > 



1 



\M\' 



4 



pp ppn 

pp pA'^ ppn^ 
pp pN^ ppn^ 



Clebsch-Gordon coefficients: 



2 
^3 



pn : — ' 



1 

^3 



^pp^p n 7r+ 



p n^ppiT 



^p n^p n TT^ 



2 1 

2^3/2 + 3^1/2 

10 2 

1 1 

2^3/2 + 3^1/2 

4 1 

3^^3/2 + 2^1/2 
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The NN NNp-cToss section 



fitting matrix elements for resonances tliat couple to the 
pp ppp^ data (similarly to l7r-cross section) 

_ 2 1 
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Resonance properties 
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The NN NNtitv-ctoss section 
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